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We calculated the Coulomb drag contribution to the resistivity and the transresistivity for a 
double layer system, in which the passive layer is a Wigner crystal pinned by impurities, and the 
active one is a metal. We found that in quasi one dimensional systems both quantities are suppressed 
at low temperatures by the constraints of energy and momentum conservation in electron-phonon 
scattering. In two dimensions both quantities are oc T**, which is consistent with Bloch law for 
the contribution of electron-phonon interaction to the resistivity. Strong impurities in the Wigner 
crystal decrease the transresistivity significantly when their number reaches the logarithm of the 
number of sites in the crystal. In contrast, the drag contribution to the resistivity is found to be 
independent of the impurities, at least as long as their density is kept zero in the thermodynamic 
limit. 



I. INTRODUCTION 

Coupled quantum wells are useful tools for studying 
electron-electron interaction in low-dimensional systems 
both experimentally and theoretically. By changing the 
distance between the layers, one can allow or prohibit 
interlayer tunneling and also control the interaction be- 
tween electrons in different layers. Moreover, by changing 
parameters such as the density of electrons, the temper- 
ature and the magnetic field one can examine different 
electronic phases. 

One effect, which has attpacted much attention in re- 
cent years, is Coulomb dragJcl, where current Ii in one 
layer ("the active layer" ) induces a voltage V2 in the other 
("passive") layer, which is kept in a zero current state 
{I2 = 0). The measured quantity is the drag resistiv- 
ity, or trans-resistivity, which for a square sample is de- 
fined as pe, = — ^ • Since electron-electron collisions have 
only indirect consequences for transport properties in a 
single quantum well, because they conserve momentum, 
the Coulomb drag is unique in directly measuring effects 
of electron-electron interaction through a transport mea- 
surement. 

Coulomb drag has been studied extensively in the 
regime in which carriers in both layers form Fermi liq- 
uids. For low enough density the two-dimensionaLelec- 
tron liquid should condense into a Wigner crysta£3. A 
clean Wigner crystal should have no electrical resistivity. 
However, it is well appreciated by now, fippm analogous 
behavior in charge density wave systemsEl, that an ar- 
bitrarily small disorder potential should pin the crystal 
so that there are no charge carriers that can flow in re- 
sponse to an arbitrarily small electric field. The disorder 
also causes the crystal to deform, so that it has a finite 
correlation length, or domain size. As a result of pinning, 
the linear response conductivity vanishes. 

In this work we study Coulomb drag in a double layer 
system in which carriers in the passive layer form a 
pinned Wigner crystal. The pinning of the crystal is due 



to its interaction with strong impurities, whose density 
is assumed to vanish in the thermodynamic limit. The 
common methods for calculating drag resistivity are noi 
directly applicable to this case: the Boltzmann equationQ 
can be used to calculate the flow of momentum from the 
active to the passive layer, but does not distinguish be- 
tween momentum flow to the electrons in the Wigner 
crystal and momentum flow to the impurities that pin 
the crystal. This distinction is important for separating 
between two different measurable quantities: the contri- 
bution of the passive layer to the resistivity of the active 
layer and the trans-resistivity, proportional to the voltage 
developing on the passive layer. The Kubo formalismErEI 
is useful for calculating the iT&ns conductivity. However, 
the inversion of the conductivity matrix, needed to ob- 
tain the transresistivity, is delicate, due to the vanishing 
zero temperature conductivity of the passive layer. Thus, 
this problem necessitates a different method. 

To that end, we introduce an harmonic potential of fre- 
quency loq, that confines the passive layer and does not 
allow it to move. When momentum flows from the active 
to the passive layer, the passive layer is displaced from 
its equilibrium position. Its displacement, which we call 
Uo, is determined by the balance between the momentum 
flux from the active layer and the restoring force exerted 
by the harmonic conflning potential. We calculate Uq, 
from which we deduce the restoring force and the tran- 
sresistivity. The introduction of the harmonic potential 
into the problem allows us to calculate the transresistiv- 
ity directly by using standard linear response formalism. 



II. THE MODEL AND THE PROBLEM 

The model is deflned as follows. There are two con- 
ductors (wires in the quasi- ID case and layers in the 
2D case). Current flows in the active one, which is in a 
weakly disordered Fermi liquid state. No current flows in 
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the passive conductor, which is in a Wigner crystal state 
pinned by impurities. 

The Hamiltonian of the system is given by 



H — Hi + H2 + 



(1) 



Here Hi and H2 are the Hamiltonians of the active and 
passive conductors, and i?i_2 is the interaction between 
them. 
Thus 

P 

where Hdisorder is the standard interaction of electrons 
in a metal with dilute impurities; the Hamiltonian of 
the passive conductor, H2, describes a crystal interacting 
with impurities: 



H2 = 



dx I p(x)^ 



So 



2m2 



' ^ iAGG2sou(x)2,5(x - R,™p) J , 

imp / 



(2) 



where u is the displacement field, p - its conjugate mo- 
mentum, So is the size of a unit cell of the Wigner crys- 
tal, 7712 - the mass of electrons constituting the Wigner 



crystal and Vs - the sound velocity of the crystal (for sim- 
plicity we assume it to be equal for the longitudinal and 
transverse modes). G is a basis vector of the reciprocal 
Wigner lattice, Aq - a parameter, defining the strength 
of the impurities in the crystal. We expect the induced 
field (which enters into the definition of trans-resistivity) 
to be independent of wo, and at the end of the calculation 
is taken to zero. In our analysis of the Hamiltonian 
(^ we assume that the impurities in the passive layer are 
dilute and strong, in a sense that is defined below. 

The HpriAratinn of this Hamiltonian is similar to one 
given inl30. We start from a discrete harmonic lattice, 
interacting with impurities: 

^2 = E + ^ E "(^OD(x, - x,)u(x,) 



V{yii + u(xj) - Rimp). 



(3) 



Here D is the dynamical matrix of the crystal, V - the 
impurity potential and the third sum extends over lattice 
and impurity sites. Now we write 



u = u 



u 



where vP gives the static deformation of the crystal due 
to the impurities, and - fluctuations around the new 
equilibrium. Expanding everything in small quantities 
to second order, we get: 
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^2 = E ^ + \ E "°(-OD(x. - x,)u°(x,) + y(x, + u«(x. 

+ ?E"'(^»)^(^»~^J>'(^J-)+ E J(u'(x^)V)(ui(xOV)V^(x, + u"(x,) 



(4) 



where there is no linear in U"'^ term, because u'^ minimizes 
the potential energy, including the impurities, and the 
crystal is harmonically pinned at each impurity site. We 
now assume that the impurity strength is much larger 
than elastic forces in the crystal, so that the lattice 
adjusts itself in order to minimize the interaction en- 
ergy with impurities at [Cach impurity site. This is a 
strong pinning conditionE2l. It is valid when impurities 
are strong and dilute. As derived in detail below, the 
impurities are strong when 



G'^Aq » 2m2VsU;o/ao 
in quasi ID, and 

2TTin2Vg 

So H—) 



G^Aa » 



(5) 



(6) 



in 2D. Here Oq is the lattice spacing. 

Assuming strong pinning, can be found simply by 



site. Assuming now that the range of the impurities is 
much larger than the lattice spacing, but much smaller 
than the inter-impurity distance, we can write 

E V"(x,+u°(x,)-R„„p) « J2 V'(x,+u"(R„„p)-R„,p), 



which is periodic in u''(Rimp) with lattice period. Ex- 
panding in a Fourier series and taking only the lowest 
Fourier components, we may approximate 



Enx,: + u°(xO-R. 

i 

~ E COs(G(u°(Rj^p) - R^mp)), 



(7) 



G 



which for attractive impurities (Ac < 0) is minimized 



when cos(G(u (R. 



■imp ) 



Rimp)) = 1- Substituting this 



minimizing V{iii + u*'(xi) — Rimp) at each impurity approximation in the Hamiltonian (H), we get: 
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^(ui(xOV)(ui(x,;)V)V^(x, + u°(x,) - R™p) 

i,imp 

« -Ug J2 (u^(x.)G)(u1(xOG)5.,_k._. (8) 



We consider both quasi ID coupled conductors (by which 
we mean wires of several conduction channels, in which 
electrons form either a Fermi liquid or a Wigner crystal), 
and two dimensional conductors. We focus mostly on the 
latter. 



We now omit the non-dynamical parts of (U) (i.e., 
the parts that depend only on u"), add a term 
5 X] "^2t^o('^(x))^ ^'^ represent a harmonic confining po- 
tential for the passive layer, and go to the continuum 
limit. We obtain the Hamiltonian (^. Since this Hamil- 
tonian does not depend on u*^, we can simplify the no- 
tation by omitting the superscript from and calling it 
u . It should be remembered, however, that from now 
on u is not the deviation of the electron position from a 
lattice site, but its deviation from the position it holds in 
equilibrium (which is shifted from the lattice site by the 
impurities). 

The interaction between the conductors is given by 



Hl-2 = Sq+G + «'5k+G-qqUk 



hG-q(quki)(qukJ ]n2Ucicl_ 



(9) 



where 712 is the density of lattice sites in the Wigner 
crystal. This is just the inter-layer Coulomb interaction 
expanded to second order in terms of the crystal dis- 
placement field. We omitted here a term arising from 
the static deformation of the Wigner crystal due to im- 
purities. Uq is the Fourier component of the screened 
interlayer Coulomb potential. We take it to be constant 
at wave length smaller than the interlayer spacing, while 
at larger wavelengths it decreases exponentially. 

The problem is posed as follows: given a current den- 
sity ji in the active conductor and j2 = in the passive 
conductor, what are the electric fields Ei + SEi and E2 in 
the conductors? (Here Ei is the field which would be in 
the active conductor in the|-absence of the passive one). 
Using the resistivity matrixQ, we can write this as 



El + SEi 
E2 



p + Sp -pD 
-PD PW 



(10) 



where p is the resistivity of the active layer in the absence 
of drag. What we need to find are 5p and po- 

We find the field E2 by calculating uq, since the two 
are related by 



E2 



7712^5^0 



(11) 



while (5Ei is found by calculating Sji— the negative con- 
tribution to the current in the active conductor at a given 
field El due to the drag. Then, for weak drag. 



6ii = -SEi/p. 



(12) 



III. RESULTS 

In quasi ID conductors the transresistivity and the 
drag contribution to the resistivity are suppressed at low 
temperatures in the limit of vanishing ujq and infinite L, 
since in that limit momentum and energy cannot be con- 
served simultaneously for electron- phonon scattering. 

In the 2D case the transresistivity is found to be 



Pd ^ Z 



mi 



7712 e^fC' V f sV^rii 
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imp 



2ti ln(L/ao) 



(13) 



where Uq is the zero wave vector component of the inter- 
layer interaction, Z - a numerical factor given by 



Jo sinh (z) 

mi is the electron mass in the active (metallic) layer, 777-2 
is the electron mass in the Wigner crystal, Vg is the sound 

velocity in the Wigner crystal, Vfs = 



j — Vg, with Vf 

being the Fermi velocity in the metallic layer (we assume 
Vf > Vs and that Vf,Vs and Vfs are of the same order of 
magnitude) ; 711 is the density of electrons, Nimp is the 
number of impurities, and L/aQ - the ratio between the 
size of the system and that of a unit cell. The impurities 
are supposed to be dilute enough, so that a condition 
Nirnp ^ ln(-Zv/ao) is satisfied. 

We see that impurities in the passive conductor de- 
crease the transresistivity. Due to the force they exert on 
the crystal, part of the momentum flux from the active 
layer is transferred to the impurities. This part does not 
lead to drag. Stretching the approximation to its limit, 
we find that the influence of the impurities becomes sig- 
nificant once their number reaches the logarithm of the 
number of electrons in the Wigner crystal. The impu- 
rity strength does not appear in the final expression as 
long as it is much larger than elastic forces in the Wigner 
crystal. The temperature dependence is consistent with 
the 2D extension of the Bloch law for the resistance due 
to electron-phonon scattering. 

The drag contribution to the resistivity in 2D is, in the 
limit of dilute impurities. 



5p 
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ni 



PD 1 



Z 



7l2777i 



27rln(L/ao) 



(14) 



711TO2 e^fC'vfsvlni 
This contribution is not suppressed by the impurity in- 
duced factor ( 1 — „ i^/r"? — T 1 . The momentum, that 

2irln(L/ao) J ' 
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flows to the impurities, does increase the resistivity in 
the active layer, although it does not induce voltage on 
the electrons in the passive layer. In the absence of im- 
purities, the ratio between the transresistivity and the 
drag contribution to the resistivity is 7ii/ri2, as required 
by the Galilean invariance. Due to this invariance the 
two electronic systems can flow at equal speeds without 
any friction resulting from their interaction. Thus when 
ji/j2 — ni/n2 we must have 

SE ^ jiSp - j2PD = 0. (15) 



where we have assumed isotropy. Then in imaginary time 
the correlator is given by 

Uu{iLu) = / Te'"^(ruo(T)p^). (18) 
Jo 

Here /3 is the inverse temperature and is the com- 
ponent of the total momentum of electrons in the active 
conductor in the direction of the applied field. This cor- 
relator is then calculated to lowest order of inter-layer 
interaction. The corresponding diagrams are shown in 
Fig. |l|. Note that the last pair of diagrams vanishes, 
since the transferred momentum is zero. 



IV. CALCULATION 

As explained in the introduction, the calculation is 
done using the linear response formalism. We first con- 
sider the trans-resistivity of the passive conductor and 
then the contribution of inter-layer interaction to to the 
resistance of the active conductor. 



A. The transresistivity 

From the Kubo formula we find the displacement of 
the Wigner crystal: 




p+G 



k+G 





!0) 



(0) 

1) ' 



p+r 



k+G 





«) = lim ^ 



d^exp^-*(K(^),j^? 



(16) 



where jo — e/m qc^Cq is the current operator for the 
active layer and E - the applied-electric field (in the same 
layer). Proceeding as usualEjO, we define the correlator 

n„ by 



{uo{uj)) = -n„(a;)— E, 
LO mi 



(17) 



n„(w^o)= J2 



kik2P 



FIG. 1. Diagrams contributing to the drag coefficient to 
second order in the electron-phonon potential. The full lines 
are electron Green's functions, the long dashed lines are 
phonon Green's functions, the short dashed ladder represents 
a difluson, and the wavy line is a current vertex. The cross 
represents a periodic external potential. 



Calculation of the topmost pair of diagrams yields an 
expression for n„: 



J^J ^?iwrpnGp(ei)p/3n^^(ei)n;^(-ei-e2)nB(e2)(i^k.k(e2)-I?kk.(f2))3Gp+q(ei+e2), (19) 

where i n the negle ct of Umklapp processes q = k = ki + k2 and ki, k2,k are in the Brillouin zone. The frequency 
LUk = \/ w^k^ + cOq is the phonon frequency for a clean crystal, S is the area of the sample, T is an electron current 
vertex function in the active layer, G is the conventional electron propagator in the active layer, and finally D is the 
phonon propagator including impurities, which should be calculated. Substituting for G, which gives delta-functions, 
integrating over ei and 62, and assuming no Umklapp (since the potential Uq is suppressed at large q), so that q = k, 
we get 



n„(w ^ 0) 



:W/3n_F(ep)nF(-ep+k) 



,4^ '° \'2m2j S'^^wo^ki^k^Wk^ 

kik2P ^ ' ^ 

xnB(ep+k - ep)Ttfip^ (£'k2k((ep+k - £p)/^) - ■Dkk2((ep+k - <^p)/^)) 



(20) 



4 



where Tt is the electron transport scattering time for the active conductor and ep - electron energy for the active 
conductor. Now we sum over p, assuming T <C and 1/d <C kf (where d is the spacing between the conductors and 
kf the Fermi wave vector). Then 



n„(c. ^ 0) - ^ -^A^(0)Z?k,o(0) 



kik2 



\2m2 



2 

hujf3 



^ti^oWki^ka^k 4 sinh^ {(ihvfk cos 6'/2) 



(^k2k(t'Ffccos0) - Dl^^ivpk COS 9)) , (21) 



where A^(0) is the extensive density of states for the 
metallic conductor. To proceed further, we have to cal- 
culate the phonon propagator D in the presence of im- 
purities. Here we calculate the phonon propagator in 
different dimensions, assuming that strength of the im- 
purities is much larger than typical clastic forces in the 
crystal, i.e., strong pinning. To formulate quantitatively 
this condition, we consider the static deformation of the 
Wigner lattice due to all impurities. In our derivation 
of H2 above, we assumed that the fragment of the lat- 
tice, that interacts with each impurity, is positioned in 
a way, that minimizes this interaction. In doing so, we 
neglected elastic forces on the scale of inter-impurity dis- 
tance and found out that, up to a lattice constant, the 
displacement of a fragment interacting with an impurity 
at Rimp satisfies vP = Rimp- The elastic forces between 
fragments will modify this displacement. We now ex- 
press this modification in terms of phonon propagator. 
In coordinate representation the propagator is 



i)(w,Xi,X2) = 



|e-*(-^0(^)([u.,(^),^.x.(O)])), (22) 



and it can be written in momentum representation as: 

EhU, 1 (', ,N„jkixi-ik2X2 
"■i^ki,k2 , , 
?777-Z — 7=== ' (23) 



ki,k2 



2A^2m2^WkiWk2 



where N2 is the number of sites in the Wigner lattice. 
The propagator D relates a point force / applied at X2 
to the displacement it generates at xi 



Mxi = -D{UJ = 0,Xi,X2)/. 



(24) 



The force exerted on the lattice by an impurity at R^ is 
(assuming G(u°(Ri) - R^) < 1) 

Consequently, the displacement satisfies the equa- 
tion 

< =-E^-^(^ = 0'^fc'^^)^'^G(R»-u^J, (25) 



where [lu — 0, xi, X2) is the propagator of a clean crys- 
tal in the coordinate representation. When the impurities 
arc dilute, —D^{uj — 0,Xfe,XA;) >> —D^{uj — 0,XA;,Xi) 
for i ^ k. Under this condition we may retain only the 



term with i = k in the sum, and the strong pinning con- 



dition, requiring that m'' >> Vi-, 



jP becomes: 



(26) 



Physically, the neglect of the i ^ k terms in the sum (|2£_ 
implies that the local displacement is determined by 
a balance between the potential exerted by the closest 
impurity and the confining harmonic potential charac- 
terized by loq. 

In ID the condition (Bq) is 



G'^Ag » 2m2VsUJo/ao, 



while in 2D: 



G^Ag » 



2TTm2v1 
soln(^)' 



(27) 



(28) 



where so is the area of the unit cell. 

The calculation of the phonon propagator in the pres- 
ence of the impurities goes as follows. We start from a 
differential equation satisfied by the propagator in the 
coordinate representation: 

{uj' + vlV^ -ojI + Y ^2^^5(x - xol D(c., X, x') 



m2 



so<5(x - x') 



m2 



(29) 



where the sum is over the impurities. This equation can 
be solved in terms of the free propagator 

l)(w,x,x') = I)^'(cj,x-x') - 

; G''AgD^\uj, X - ^,)b{uj, x„ x'), (30) 



which gives 

^(a;,xfe,x') = V(^"')fc,^^(cc',x, -x'), (31) 
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with matrix A defined by 

A,, = + G^AcD^iiJ, X, - x,). (32) 
Substituting this back into (pO|), we obtain 
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L»(w,x,x') = (w,x,x'; 



E 



Using 



we can obtain an expression for the propagator in the momentum representation: 



h 



— 'tk;iXi+'tk2Xj 



SWkifJkika 



E 



g-ikiXigjk2Xj 



hrriyNy 



(33) 



(34) 



where the free propagator in the momentum representation was used: 



^k.M^-^k.k.^la 



2cjki 



M(c. + z<5)2-0- 

In the Hmit of strong pinning, the propagator satisfies an important relationship: 



k 



ikixi 



Dk^ojto ^ 0) 

V^oWki 



1 



G^AcDfiuj = 0,x = 0) 
I 



(35) 



(36) 



where x; is an impurity site. This relationship means 
that a force, when applied at an impurity site, causes a 
much smaller deformation, than when applied anywhere 
else, since it is opposed by the impurity. The proof of 
this relationship follows from the equation satisfied by 
the propagator D, namely D = A~'^D^. In the strong 
pinning limit Aij ~ G^AGD^{u!,x.j — x^) and hence (36). 



From the above relationship it follows that for all 



ki 



(^kk.(^^))* 



0, 



(37) 



where = D — Df is a. part of the propagator due to 
the impurities. To see this, we note using Eq. (34) that 
-Dkio(0)/^Wki is a sharply peaked function of ki around 
ki = 0: its free part is proportional to (5ki.0i ^^nd the 
part due to impurities is peaked around zero with width 
uq/vsi because of the denominator l/w^^^. In contrast, 
^k2k(^) varies slowly in that region, except for the ex- 
ponential factor e~'^^^^^ = g-^'^^i g^'^i^i ^ Hence the sum- 
mation is actually of the form (p6[). 



Using the last relationship in (pT|), and taking into ac- 
count that the typical energy transfer in an inter-layer 
scattering event, vpkcosO, is of the order of the tem- 
perature, and thus much larger than wq, we see that 
the effect of impurities enters Ec^. ( |2l| ) only in the first 
phonon propagator, where the frequency is zero, while 
for the second phonon propagator, that carries a non- 
zero frequency, only the free part of D contributes. Thus, 



-iV(0)^i?«„(0)( — )2- 



xhuj/STtlq. 



Vfhk cos 9 



Asmh-'{f3vfhkcoa0/2) 



: (X]T%^'5(w/fc cos 61 TWfe)). 



(38) 



To proceed further with the calculation, we should con- 
sider different dimensions. In ID (where electrons in the 
active wire are in the Fermi liquid regime) cos 9 = ±1, 
so that the momentum conservation condition is never 
satisfied (in the limit of vanishing luq and vanishing tem- 
perature) . This prevents the crystal from getting any mo- 
mentum from the conductor and thus the trans-resistivity 
is zero in the low temperature limit. 

The situation is different in two dimensions. In order 
to calculate n^, in 2D we have first to find I?oo(0). For 
this we write approximately eq.(B2h for a; = 0: 



A 



j,iuj = 0) = G'Ac 



i?^(cj = 0,xj-x,). 



(39) 



Now we assume that the number of impurities is so low 
that the distance between them is of the order of the 
system size (hereby neglecting fluctuations in which im- 
purities could be close together). Then the matrix A 
satisfies 



Aijiuj = 0) « Auiuj = 0) 



(40) 



for i ^ j. The matrix A then can be inverted approxi- 
mately: 



1 



G^AaiDfiu = 0,x = 0))2 
xL)-^(c^ = 0,Xj-x,), (41) 



6 



where S{i^) = 1 for i = j and —1 for i ^ j. Substituting 
this into Q, we find I^oolO) in 2D: 



(42) 



where Nimp is the number of impurities. The ln(^;^^) 

originates from the expression for D fiuj — Q,yii = 0, X2 = 
0) in 2D. 

Substituting expression ( ^2| ) for the propagator into 
Eq.(p8[) and choosing 



we obtain the transresistivity Eq. (|l3|). 

As evident in Eq. ([l3|), the impurities decrease the 
trans-resistivity as they absorb part of the momentum 
transferred from the active layer. That part of the trans- 
ferred momentum does not generate a drag voltage, and 
therefore does not lead to trans-resistivity. Under the 
strong pinning approximation, the influence of the impu- 
rities becomes significant when their number is compara- 
ble to ln(L/ao), even though their density is still zero in 
the thermodynamic limit. 



B. Drag contribution to resistivity 

We now turn to the way electron-electron scattering 
between different layers affects the resistivity of the ac- 
tive layer. We again apply the Kubo formula to find, to 
second order in screened inter-layer interaction, the cur- 
rent density in the active layer as a response to an electric 



field in that same layer (we expect it to be negative as 
the drag should decrease the current): 



fiuoS 



dte^-\[m,f,])^f^. (43) 



Introducing the correlator from 

j" = - lim 3nf (w)E'3, 
we have in imaginary time: 



n(icj) = - 



dre--(Tjg(T)jg(0)), 



(44) 



(45) 



where we again assumed isotropy. The corresponding di- 
agrams are shown in Fig. 0. 






FIG. 2. Diagrams contributing to reduction of conductivity 
of the active conductor 

Calculating them, we obtain 



n = -(^)^E^^ 

mi ^ \ S 

p,k 



"-2 



"■f (£p) - »F(-£p+k) 
2m2Wfe~"4sinh^(/3(ep+fc - ep)/2) 



2u2 



(ep+fe - ep)nr/k: 



(46) 



In one dimension, when the passive layer is clean, the 
constraints of energy and momentum conservation to- 
gether with the Pauli principle suppress the resistivity 
at low temperature. In the presence of impurities in the 
Wigner crystal momentum conservation is relaxed, and 
we obtain a contribution, which is inversely proportional 
to the length of the system. This contribution can be 
neglected, since in the macroscopic limit we assume the 
length to go to infinity, while the density of impurities is 
kept zero. 

In two dimensions, in the absence of impurities mo- 
mentum and energy can be conserved simultaneously, 
and inter-layer interaction yields the following contribu- 
tion to the current 



h^mim2VfsV^ 



ZEi = (5ctEi, 



(47) 



from which the drag contribution to the resistivity is eas- 
ily found: 



Sp = —p^Sa = 



(48) 



In the presence of impurities the phonon propagator is 
composed of a free part and an impurity- induced part. 
The contribution of the former remains Eq. (p8|), while 
the contribution of the latter is again inversely propor- 
tional to the system size. Hence in the macroscopic limit 
the impurities do not affect Eq. (^8|). 

When the passive layer is clean, the ratio between the 
trans-resistivity and the drag contribution to resistivity 
of the active layer is ni/n2. This is consistent with the 
Galilean invariance requirement. 
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V. SUMMARY 



We considered a bi-layer system in which the passive 

layer is a Wigncr crystal pinned by impurities and the 
active layer is a conductor, and calculated the transre- 
sistivity between two layers and the drag contribution to 
the resistivity of the active layer. We focused on the case 
of a small number of strong impurities, whose density 
vanishes in the thermodynamic limit. We found that in 
quasi ID case both quantities are exponentially small at 
low temperature due to constraints of energy and mo- 
mentum conservation. 

In two dimensions we found a T^-dependence of both 
quantities, consistent with the Bloch law for the electron- 
phonon interaction contribution to resistivity. Impurities 
in the Wigner crystal decrease the transresistivity signif- 
icantly when their number reaches the logarithm of the 
number of sites in the crystal. In the limit of strong pin- 
ning the transresistivity and the active layer resistivity 
are independent of the impurity strength. 

While a minute density of impurities is sufficient to 
affect the transresistivity, we find that it does not signif- 
icantly affect the contribution of the Wigner crystal to 
the resistivity of the active layer. 
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